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Abstract 

We present a measurement of the Standard Model CP violation parameter 
sin 2^1 based on a 29.1 fb _ data sample collected at the Y(4£>) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e~ collider. One 
neutral B meson is fully reconstructed as a J/ipKs, ip(2S)Ks, XciKs, r/ c Ks, 
J /i/)Kl or J/tpK* decay and the flavor of the accompanying B meson is 
identified from its decay products. From the asymmetry in the distribution 
of the time intervals between the two B meson decay points, we determine 
sin 20i = 0.99 ± 0.14(stat) ± 0.06(syst). We conclude that we have observed 
CP violation in the neutral B meson system. 
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Kobayashi and Maskawa (KM) proposed, in 1973, a model where CP violation is incor- 
porated as an irreducible complex phase in the weak- interaction quark mixing matrix [|TJ. 
The idea, which was presented at a time when only the u, d and s quarks were known to 
exist, was remarkable because it required the existence of six quarks. The subsequent dis- 
coveries of the c, h and t quarks, and the compatibility of the model with the CP violation 
observed in the neutral K meson system led to the incorporation of the KM mechanism into 
the Standard Model, even though it had not been conclusively tested experimentally. 

In 1981, Sanda, Bigi and Carter pointed out that the KM model predicted large CP 
violation in certain decays of B mesons for a range of quark mixing parameters. Subsequent 
measurements of the B meson lifetime || and the discovery of B°B° mixing || indicated 
that the parameters lie within such a range. Thus, measurements of CP violation in B 
meson decays provide important tests of the KM model. 

The model predicts a CP violating asymmetry in the time-dependent rates for initial B° 
and B° decays to a common CP eigenstate, fcp 0. In the case where fcp = (cc)K°, the 
asymmetry is given by 



A ® = ^CtllTnoZt:] = -Osm20 lS inAm d t, 



where T(B° (B°) — > fcp) is the decay rate for B°(B°) to fcp at a proper time t after 
production, £f is the CP-eigenvalue of fcp, Am d is the mass difference between the two B° 
mass eigenstates, and 0i is one of the three internal angles of the Unitarity Triangle, defined 

as 0i = 7r — arg I _y»y d ) [§]. For the (cc)K° decays, both the ambiguity due to strong 

interactions and the contribution from direct CP violation are expected to be small [H. 

Our previous determination, using a data sample taken in 1999-2000, found sin 20! = 
0.58io;il(stat) to'io(syst) ||, which is consistent with the KM model constraints from indirect 
measurements ]7). Although the combination of this result with other measurements of 
sin 20i [[I strongly indicates violation of CP symmetry in B meson decays, the published 
results are still not conclusive. In this Letter we report a new measurement of sin 20i that 
uses improved reconstruction algorithms and incorporates data taken in 2001 to achieve a 
four-fold increase in the size of the event sample. The result reported here includes the 
earlier data and supersedes the previous value. All data samples have been analyzed and 
reconstructed with the same consistent procedure. 

We use a 29.1 fb~ data sample, which contains 31.3 million BB pairs, collected with 
the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider 0]. KEKB 
operates at the T(4S) resonance (y/s = 10.58 GeV) with a peak luminosity that exceeds 
4 x 10 33 cm _2 s~ 1 . The Belle detector is a large-solid-angle magnetic spectrometer that 
consists of a three-layer silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), a mosaic of aerogel threshold Cerenkov counters (ACC), time-of-flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside of 
the coil is instrumented to detect Kl mesons and to identify muons (KLM). The detector 



is described in detail elsewhere [ 10 1 



We measure sin 20! using B°B° meson pairs produced at the T(45) resonance. The two 
mesons remain in a coherent p-wave state until one of them decays. The decay of one of the 
B mesons at time t tag to a final state, / tag , which distinguishes between B° and B°, projects 



the accompanying B meson onto the opposite 6-flavor at t tag ; this meson decays to fcp at 
time tcp- CP violation manifests itself as an asymmetry A(At), where At is the proper 
time interval At = tcp — tt ag . At KEKB, the T(4S) is produced with a Lorentz boost of 
/?7 = 0.425 nearly along the electron beamline (z). Since the B° and B° mesons are nearly at 
rest in the T(4S) center of mass system (cms), At can be determined from the displacement 
in z between the fcp and / tag decay vertices — i.e. At ~ (z C p — z ta:g )/ [3^c = Az/fljc. 

The measurement requires the reconstruction of B° — > fcp decays, the determination of 
the 6-flavor of the accompanying (tagging) B meson, the measurement of At, and a fit of 
the expected At distribution to the measured distribution using a likelihood method. 

We reconstruct B° decays to the following CP eigenstates |TT| : J/i/jKs, if)(2S)K$, XaKs, 
rj c Ks for £/ = — 1 and J/ipKi for £/ = +1. We also use B° — >■ J/ipK*° decays where 
K*° — > Ksii . Here the final state is a mixture of even and odd CP, depending on the relative 
orbital angular momentum of the J /if) and K*°. The CP content is determined from a fit to 
the full angular distribution of all J/if)K* decay modes other than K*° — > Ksn . We find that 
the final state is primarily £) = +1; the £/ = -1 fraction is 0.19±0.04(stat)±0.04(syst) pf. 

J /if) and if>(2S) mesons are reconstructed via their decays to £ + £~ (£ = //, e). The if)(2S) 
is also reconstructed via J/ip7i + 7i~, and the Xci y i a J/i 3 !- The r\ c is detected in the K + K~7r° 
and KsK~n + modes. For the J/if)Ks mode, we use K$ — » 7r + 7r~ and 7r°7r° decays; for other 
modes we only use Ks — > 7r + 7r~ . 

The J/ip, if>{2S) and Ks selection has been described elsewhere 0. For XaKs decays, 
we select \c\ -^ J/ipJ decays, rejecting 7's that are consistent with 7r° — >■ 77 decays, and 
use the requirement 385 < M 7 « — M« < 430.5 MeV/c 2 . For r\ c decays, we distinguish kaons 
from pions using a combination of CDC energy loss measurements, flight times measured in 
the TOF, and the response of the ACC. Candidate r\ c -^ K + K~n° (KsK~it + ) decays are 
selected with a KKtt mass requirement that takes into account the natural width of the r\ c . 
For J/ipK* (Ks7T°) decays, we use Ksn combinations that have an invariant mass within 
75 MeV/c 2 of the nominal K* mass. We reduce background from low- momentum 7r°'s by 
requiring cos 9k* < 0.8, where 9k* is the angle between the Ks momentum vector and the 
K*° flight direction calculated in the K*° rest frame. 

We identify B decays using the energy difference AE = E c ^ s — E^ m and the beam- 



energy constrained mass Mb c = y(E^ m ) 2 — (p^ 1 * 3 ) 2 , where E^ m is the cms beam energy, 
and E^ 18 and p c ™ are the cms energy and momentum of the B candidate. 

Figure [l] shows the combined M^ c distribution for all channels other than J /ipK^ after 
a mode- dependent requirement on AE. The B meson signal region is defined as 5.270 < 
M bc < 5.290 GeV/c 2 . Table | lists the numbers of observed candidates (iV ev ) and the 
background (iVbkgd) determined by extrapolating the rate in the non-signal AE vs. M^ c 
region into the signal region. 

Candidate B° -^ J/ipKi decays are selected by requiring ECL and/or KLM hit pat- 
terns that are consistent with the presence of a shower induced by a neutral hadron. The 
centroid of the shower is required to be in a 45° cone centered on the Kl direction that 
is inferred from two-body decay kinematics and the measured four-momentum of the J/if). 
We reduce the background by means of a likelihood ratio that depends on the J /if) cms 
momentum, the angle between the Kl and its nearest-neighbor charged track, the charged 
track multiplicity of the event, the extent to which the event is consistent with a B + — > 
J/if)K* + (KL7r + ) hypothesis, and the polar angle with respect to the z direction of the re- 



constructed B° meson in the cms. In addition, we remove events that are reconstructed as 
B° -> J/ipK s , J/ipK*°(K + 7r-,K S 7r°), B + -» J/^ + , or J/iPK* + (K + tt° } K s ir + ) decays. 
Finally, Kl clusters with positions that match photons from reconstructed 7r°'s are also 
rejected. 

Figure [^ shows the p c ™ distribution, calculated with the B° — > J/ipK L two-body decay 
hypothesis. The histograms are the results of a fit to the signal and background distributions. 
The shapes are derived from Monte Carlo simulations (MC) |13| , and the normalization and 
peak position of the signal are allowed to vary. There are 397 entries in the 0.2 < p c ™ < 
0.45 GeV/c signal region with KLM clusters. There are 172 entries in the range 0.2 < 
p c g is < 0.40 GeV/c with clusters in the ECL only. The fit finds a total of 346 ± 29 J/^K L 
signal events, and a signal purity of 61%. 

Leptons, charged pions, and kaons that are not associated with a reconstructed CP 
eigenstate decay are used to identify the flavor of the accompanying B meson. Initially, the 
6-flavor determination is performed at the track level. Several categories of well measured 
tracks that distinguish the 6-flavor by the track's charge are selected: high momentum 
leptons from b — > c£~V, lower momentum leptons from c — > sl + v, charged kaons and A 
baryons from b —>■ c — > s, high momentum pions that originate from decays of the type 
B° — > £)(*)- (7r + ,p + , a+,etc.), and slow pions from D*~ — > D°ir~. We use the MC to 
determine a category-dependent variable that indicates whether a track originates from a 
B° or B°. The values of this variable range from —1 for a reliably identified B° to +1 
for a reliably identified B° and depend on the tagging particle's charge, cms momentum, 
polar angle and particle-identification probability, as well as other kinematic and event shape 
quantities. The results from the separate track categories are then combined to take into 
account correlations in the case of multiple track-level tags. This stage determines two 
event-level parameters, q and r. The first, q, has the discrete values q = +1 when the 
tag-side B meson is more likely to be a B° and —1 when it is more likely to be a B°. The 
parameter r is an event-by-event flavor-tagging dilution factor which ranges from r = for 
no flavor discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort 
data into six intervals of r, according to flavor purity; the wrong-tag probabilities for the 
final fit are determined from data. 

The probabilities of an incorrect flavor assignment, Wi (I = 1,6), are determined di- 
rectly from the data for the six r intervals using exclusively reconstructed, self-tagged 
B° -> D*-£+v, DW-tt + , D*-p+ and J/ipK*°(K+ir-) decays. The 6-flavor of the accompa- 
nying B meson is assigned according to the flavor-tagging algorithm described above. The 
exclusive decay and tag vertices are reconstructed using the same vertexing algorithm that 
is used in the analysis to measure CP asymmetry. The values of wi are obtained from the 
amplitudes of the time-dependent B°B° mixing oscillations: (iVoF — Nsf)/{Nof + Nsf) = 
(1 — 2wi) cos(Am,dAt). Here iVoF and iVsF are the numbers of opposite and same flavor 
events. We fix Arrid at the world average value |I3 |. Table [Tl] lists the resulting wi values 



together with the fraction of the events (fi) in each r interval. The total effective tagging 
efficiency is £, /,(1 - 2w t ) 2 = 0.270 ± 0.008(stat)+J5;JJg(syst). 

The vertex positions for the fcp and / tag decays are reconstructed using tracks that have 
at least one three-dimensional coordinate determined from associated r-<p and z hits in the 
same SVD layer along with one or more additional z hits in the other layers. Each vertex 
position is required to be consistent with the interaction point profile smeared in the r-0 



plane by the B meson decay length. The fcp vertex is determined using lepton tracks from 
J /if} or ip(2S) decays, or prompt tracks from r\ c decays. The / tag vertex is determined from 
well reconstructed tracks not assigned to fcp- Tracks that form a Ks are not used. The 
MC indicates that the typical vertex-finding efficiency and vertex resolution (rms) for zcp 
(^tag) are 92 (91)% and 75 (140) jum, respectively. 

The proper-time interval resolution for the signal, R S i g (At), is obtained by convolving a 
sum of two Gaussians (a main component due to the SVD vertex resolution and charmed 
meson lifetimes, plus a tail component caused by poorly reconstructed tracks) with a function 
that takes into account the cms motion of the B mesons. The fraction in the main Gaussian 
is determined to be 0.97 ± 0.02 from a study of B° -> D*-tt + , D*~p + , D-tt + , J/ipK*°, 
J/if>K s and B + — > D°n + , J/ipK + events. The means (// ma i n , /-ttaii) and widths (cr ma _ in , 
°"taii) °f the Gaussians are calculated event-by-event from the fcp and / tag vertex fit error 
matrices and the \ 2 values of the fit; typical values are /i ma i n = —0.24 ps, /x ta ii = 0.18 ps 
and (T ma i n = 1.49 ps, Otaii = 3.85 ps. The background resolution Rbk a (At) has the same 
functional form but the parameters are obtained from a sideband region in M^ c and AE. 
We obtain lifetimes for the neutral and charged B mesons using the same procedure; the 
results |15| agree well with the world average values. 



After vertexing we find 560 events with q = +1 flavor tags and 577 events with q — — 1. 
Figure |3] shows the observed At distributions for the g£/ = +1 (solid points) and q£f = — 1 
(open points) event samples. There is a clear asymmetry between the two distributions; this 
demonstrates that CP symmetry is violated. 

We determine sin 2(f>i by performing an unbinned maximum-likelihood fit of a CP vio- 
lating probability density function (pdf) to the observed At distributions. For modes other 
than J/i/jK*° the pdf expected for the signal is 

V Big (At,q,Wi,^ f ) 

= e ~^Z B ° { 1 ~ ^f q{l ~ 2wi) Sin 2(t>1 M^rn d At)}, 

where we fix r^o and Am,; at their world average values [0. The pdf used for the background 
distribution is P b k g (At) = / r e~' A ''/ Tbk s/2r bkg + (1 — f T )5(At), where f T is the fraction of the 
background component with an effective lifetime Tbkg and S is the Dirac delta function. For 
all fcp modes other than J/i/iKl, a study using events in background-dominated regions of 
AE vs. Mbc shows that f T is negligibly small. For these modes, Pbk g (At) = 5 (At). 

The J/ipK^ background is dominated by B — > J/ifiX decays where some final states are 
CP eigenstates. We estimate the fractions of the background components with and without 
a true K L cluster by fitting the p c ™ distribution to the expected shapes determined from 
the MC. We also use the MC to determine the fraction of events with definite CP content 
within each component. 

The result is a background that is 71% non-CP modes with Tbk g = Tb- For the CP- 
mode backgrounds we use the signal pdf given above with the appropriate £y values. For 
J/ifiK*(K L 7i°), which is 13% of the background, we use the £y = —1 content determined 
from the full J/ipK* sample. The remaining backgrounds are £/ = —1 states (10%) including 
J/ipKg, and £/ = +1 states (5%) including i/j(2S)Kl, Xc\Kl and J/ipir . 

For the J /ipK* mode, we include the At and transversity angle 6 tr [|Tj| distributions in 
the likelihood |L^]. We use the £/• content determined from the full angular analysis. 



Each pdf is convolved with the appropriate R(At) to determine the likelihood value for 
each event as a function of sin 20i : 



J{f sig P sig {At', q, w h Z f )R sig {At - At') 
+ (1 - f sig )V hkg {At')R hkg (At - At')}dAt', 



where / s ; g is the probability that the event is signal, calculated as a function of p c ™ for 
J/i/jK~l and of AE and Mb c for other modes. The only free parameter is sin20i, which is 
determined by maximizing the likelihood function L = Y[ i P i , where the product is over all 
events. 

The result of the fit is 

sin 20! = 0.99 ± 0.14(stat) ± 0.06(syst). 

In Fig. ||(a) we show the asymmetries for the combined data sample that are obtained by 
applying the fit to the events in each At bin separately. The smooth curve is the result of 
the global unbinned fit. Figures |](b) and (c) show the corresponding asymmetries for the 
(cc)Ks (£/ = —1) and the J /ipK L (£j = +1) modes separately. The observed asymmetries 
for the different CP states are opposite, as expected. The curves are the results of unbinned 
fits applied separately to the two samples; the resultant sin 20! values are 0.84 ± 0.17(stat) 
and 1.31 ± 0.23(stat), respectively. 

The systematic error is dominated by uncertainties due to effects of the tails of the vertex 
distributions, which contribute 0.04. Other significant contributions come from uncertainties 
(a) in wi (0.03); (b) in the resolution function parameters (0.02); and (c) in the J/ipKi 
background fraction (0.02). The errors introduced by uncertainties in Am^ and t b o are 0.01 
or less. 

We performed a number of checks on the measurement. Table [Til] lists the results obtained 
by applying the same analysis to various subsamples. All values are statistically consistent 
with each other. The result is unchanged if we use the wi's determined separately for 
/tag — B° and B°. Fitting to the non-CP eigenstate self-tagged modes B° —> D^*^tt + , 
D*~ p + , J/tpK*°(K + TT~) and _D*~£ + z/, where no asymmetry is expected, yields 0.05 ± 0.04. 
The asymmetry distribution for this control sample is shown in Fig. f|(d). As a further 
check, we used three independent CP fitting programs and two different algorithms for the 
/tag vertexing and found no discrepancy. 

We conclude that there is large CP violation in the neutral B meson system. A zero 
value for sin 20! is ruled out at a level greater than 6a. Our result is consistent with the 
higher range of values allowed by the constraints of the KM model as well as with our 
previous measurement. 
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TABLE I. The numbers of observed events (iV e - 
signal region for each fcp mode. 



and the estimated background (iVbkgd) in the 



Mode 



N c 



IVbkgd 



J/ip(£ + l-)K s (ir+TT-) 
J/^{£ + £-)Ks{n 7r ) 
4){2S)(£ + l-)K s {-K + K-) 
ip(2S){J/ip7r + Tr~)Ks{n + 7r-) 

r] c (K + K-ir )K s (TT + TT-) 

J/i)K* Q {K s ^) 



457 
76 
39 
46 
24 
23 
41 
41 



11.9 

9.4 
1.2 
2.1 
2.4 
11.3 
13.6 
6.7 



Sub-total 



747 



58.6 



J/^{£ + £-)K L 



569 



223 



TABLE II. The event fractions (//) and incorrect flavor assignment probabilities (wi) for each 
r interval. The errors include both statistical and systematic uncertainties. 



/ 



ft 



W[ 



T465 Tlnmj 
u.4oo_ 0009 

n orro+0.015 

u.ooz_ 0014 

n 24 o+0.021 

n i 7fi +0.022 
U- 1 '"-0.017 
n11n +0.022 

u.iiu_ 0014 
041+ - 011 

u - u ^ i -0.010 



1 

2 
3 

4 
5 
6 



0.000 - 0.250 
0.250 - 0.500 
0.500 - 0.625 
0.625 - 0.750 
0.750 - 0.875 
0.875 - 1.000 



0.405 
0.149 
0.081 
0.099 
0.123 
0.140 



TABLE III. The values of sin2c/>i for various subsamples (statistical errors only). 



Sample 



sin 2(f>\ 



/tag 
/tag 



B° (q 

W (q 



+ 1) 
-1) 



0.84 ±0.21 
1.11 ±0.17 



(cc)Ks except J/iPKs(tt + tt~ 
J/tpK* (K s iT ) 



0.81 ±0.20 
1.00 ±0.40 
1.31 ±0.23 
0.85 ± 1.45 



All 



0.99 ±0.14 
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FIG. 1. The beam-energy constrained mass distribution for all decay modes combined other 
than JjijjKi,. The shaded area is the estimated background. The signal region is the range 
5.27-5.29 GeV/c 2 . 
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FIG. 2. The p§" s distribution for B° — > J/v/jKl candidates with the results of the fit. The solid 
line is the signal plus background; the shaded area is background only. The signal region for KLM 
(ECL-only) clusters is 0.2 < p c ^ s < 0.45(0.40) GeV/c. 
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FIG. 3. At distributions for the events with q£f = +1 (solid points) and q£f = —1 (open 
points). The results of the global fit (with sin2</>i = 0.99) are shown as solid and dashed curves, 
respectively. 
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FIG. 4. (a) The asymmetry obtained from separate fits to each At bin for the full data sample; 
the curve is the result of the global fit. The corresponding plots for the (b) {cc)Ks (£/ = —1), (c) 
J '/tpKi, (j£f = +1), and (d) B° control samples are also shown. The curves are the results of the 
fit applied separately to the individual data samples. 
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